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Abstract 

A functional method for calculating averages of the time-ordered exponential of 
a continuous isotropic random N x N matrix process is presented. The process is 
not assumed to be Gaussian. In particular, the Lyapunov exponents and higher 
correlation functions of the T-exponent are derived from the statistical properties of 
the process. 

The approach may be of use in a wide range of physical problems. For example, 
in theory of turbulence the account of non-gaussian statistics is very important since 
the non-Gaussian behavior is responsible for the time asymmetry of the energy flow. 
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1 Introduction 

Sets of linear stochastic differential equations appear in different physical problems related 
to quantum mechanics and held theory, turbulence, low temperature physics etc. Their 
*E-mails: asil72@mail.ru, sirota@lpi.ru, zybin@lpi.ru 
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formal solution is given by the time-ordered exponential, but calculation of the statistical 
moments is still a challenge. 

The long-time evolution of the time-ordered product of arbitrary (not Gaussian) random 
matrices in the discrete case was investigated in n. PI. and |3]. The existence of the 
Lyapunov spectrum was proved, though no recipe to calculate the Lyapunov exponents or 
other averages was given. A functional integration method to calculate the averages in the 
case of Gaussian 5-correlated processes was introduced in [1] and |5]. 

However, in many applications one needs to calculate the T-exponentials of non- 
Gaussian processes. These are, e.g., held theories with interactions. Also, all the processes 
that produce non-zero correlators of odd orders, or those in which the Lyapunov spectrum 
is not even, are knowingly non-Gaussian. This is just what occurs in the theory of magnetic 
dynamo (see 0 ) and in the theory of turbulence: the third-order correlator is related to the 
energy dissipation rate (|Z118]), and the asymmetry of the Lyapunov spectrum is necessary 
to provide the time anisotropy of a turbulent how ([9l[l0]). Besides, the velocity gradient 
tensor, which is an important object in the problems related to passive scalar (Hues]) 
and in derivation of scaling exponents (Biia), is proved to have non-Gaussian distribution 
f|14]L In this paper we propose a functional integration method for calculation the averages 
of continuous products of A^ x A^ random matrices A{t). The probability distribution of the 
matrices is assumed to be isotropic, Gaussianity is not required. The method simplihes the 
functional integrals signihcantly, in the case of ^-correlation it allows to express all the sta¬ 
tistical characteristics of the T-exponential in terms of the moments of A{t). The main idea 
of the approach is to decompose the matrix product into the rotational and deformational 
components, and then change the variables. Although the change of variables is non-local 
(relative to t), it allows to minimize the non-locality and exclude the time ordering from 
the path integral. One more important point is that we use a cumulant function to analyze 
the 5-processes. This allows to avoid renormalizations, and to get a simple expressions for 
all the averages. 
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2 T-exponential of a random process 


Let A{t), 0 < t < T be a random process taking on a valne of iV x real matrices. Its 
statistics is defined by the measnre 

N 

DAP[A\= J] l[dAkpP[A], 

0<t<T k,p=l 

where P[A] is the probability density functional. 

Consider the random matrices Q{t) that satisfy the equation 

dtQ = QA, Q (0) = i (1) 

The formal solution to the equation can be written in terms of the anti-chronological ex¬ 
ponential ([15]): 

Q{t)=Texpij A(r)dr| j dTi...dTnT {A{ti) ...A{Tn)) (2) 

Vo / ”■ 0 

where 

T {A (ri) ...A (Tn)) = A (Ti J ...A (nj , n, < ... < 

is the antichronological product operator. The alternative way to describe the solution is 
the Volterra multiplicative integral ([16]): 

t 

Q{t) = Y[{^ +A{T)dT) (3) 

T=0 

Thus, the T-exponent is equivalent to the inhnite matrix product. 

We are interested in the averages 

{Fm= J DAP[A]FIQ], (4) 

where F[Q] is some functional. These expressions contain the T-exponent which is not an 
easy object to deal with. Below we simplify them and reduce, in the case of isotropically 
distributed processes, to the path integrals of some exponents without any time-ordering. 

3 Change of variables 

We make the Iwasawa decomposition of the matrix Q: 

Q = zdR (5) 
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where z is an upper triangular matrix with diagonal elements equal to 1, d is a diagonal 
matrix, R is an orthogonal matrix: RjiRj^ = Sik- 

We recall that Q as well as z, d,R is a. function of time. From ([T]) it follows A = Q~^dtQ; 
thus, Eq. (E]) produces the decomposition of A: 

A = R^XR, X = p + C + e (6) 

where 

p = d~^dtd , C = d~^z~^ {dtz) d (7) 

e = (dtR) R^ (8) 

We see that p is diagonal, C, is an upper triangular matrix with zeroes in the main diagonal, 
and 9 is antisymmetric: 


p = diag{pi,..., Pn) , % = -9ji , Qj = 0 if i>i 


Now we consider X{t) as independent functional variables; then ([6]) should be under¬ 
stood as 

A = R^[X] X R [X] (9) 


Where R[X] is determined by ([8]): 


m 


Texp / e{T)dT 


( 10 ) 


\o / 

Note that R depends only on the ^-component of X; this will simplify further calcula¬ 
tions signihcantly. This also makes it possible to separate the rotational part of Q from 
deformation. 

The averages (jl]) can be rewritten as 


{F IQ]) = j DX J|A'] P [fl^lA'I A fl[A]] F[Q[X]] (11) 

The Jacobian J[X] will be calculated in the next section. 

The advantage of this change of variables is in dealing with a simple ’rotational’ T- 
exponent onii. instead of the complicated T-exponent (E]). As we will see below, this one 
does also vanish in the isotropic case. 
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4 The Jacobian 


In this section, we calculate the functional Jacobian 

6Aij {t) 


J = Det 


(f) 


First, to simplify the notations, we introduce the multiindices (hereafter denoted by Greek 
letters): 

ol = {ii — Aj^j , ~ Xkp ) Q:,/3,... l..A^ 

The transformation ([H]) can be presented as 


A^ = ?li^f,[X]Xp , 

where 3f?o,^[X] is the x matrix that satisfies the relation: 


( 12 ) 


3?, 


(ij),{kp) RkiRpj 


(13) 


(This corresponds to multiplying X by 3? from right and by from left). The functional 
derivative of (IT^ is equal to 


“ - 3fJ„(3 (t) S (t t') + j X^ (t) 


JX, 


(14) 


Now, as it follows from flT^ and flTOj) . 3f?(f) is determined completely by the function 6[t) 
at r < f. Thus, 3fjQ^(f) does not depend on X^{t') 

According to ([HD, the Jacobian is a block-triangular matrix: 


Consequently , its determinant is a continuous product of the determinants of the diagonal 
{t = t') X2 X X2 blocks: 


J = Det 


JAo (f) 

jx^ (to 


J {t) , J{t) = det 


0<t<T 


6Aa{t) 

SXfsit) 


(15) 


Now, let us calculate J(t). From (ITTD it follows 

/ J3f? (f) 

J it) = det (f) J (0) + {t) 
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Note that is an orthogonal matrix, since 


RkiRpj RmiRnj ^km^pn ^(kp),{mn) ^aj3 


So, det 3? = 1, and one can multiply the bracketed expression by 3?^. Then, assuming 
5 (0) = {dt)~^ and omitting the insufficient normalization multipliers, we obtain 

J {t) = det ^ 6 ai 3 + = exp (ti dt 


5X - ^ 


From firsj) it follows 


J = exp \ f G dt\ , G = 3? 


53fJ,, 


M7 Y 


(16) 


To calculate G, we return from multiindices to usual matrix notations: 

d {^RniRmj) 


G — RhiR 


ki^pj' 


_ Y 

r -^nm 

OJikp 


( 17 ) 


Now we need to calculate the variation derivative at coinciding time. The R(t) 

dependence of X(t) is determined by (E]), (ITOil . Thus, of all the components of the 
matrix X, R depends on N{N — l)/2 independent components of only. 

Since Xij coincides with 9ij as i > j, we get 

d 0 ij{t) Sikdjp if k>p 

dXkpit) 1 0 if k < p 

Hence, = 0 ii k < p. 

For the rest, we use the Volterra presentation ([3]) of the T-exponent ffHJ]) : 

0 


6 Rij{t) 


d 0 kp{t) 


k>p 


d9kp{t) 


JJ(1 + 0(r)dr) 


\T=t 


JJ(1 + 0(r)dr) 


<T=t 


dOmnir') 

d 0 kp{t) 


dr' ( (1 + 0 {T)dT) 


nj 


Further, 


d9mn{r') 


d(^T f) {dmkdnp d^ipdnk') ; 


d9kp{t) 

and we make use of the time-ordering and note that the derivative can only make a non-zero 
contribution when acting on the hrst (r' = t) multiplier. Then the rest of the multipliers 
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make R again, and 5{t — r')dT' = \ because t is the boundary of the integral. (Effectively, 
only one half of the delta-function is integrated.) The integral becomes more compact: 

0 

’ n _L 1 — _A. A .w. 

^mj 


5Rij{t) 


S0kp{t) 

Thus, 


k>p 


^ ^ip^mk^ f (1 H“ | i^^ik^mp ^ip^mk^ Rr 


\T=t 


mj 


5Rij{t^ \ 2 i^^ik^mp ^ip^mk^ Rmj i k ^ p 


dXkpit) I 0 , k < p 

As we substitute this into flT7|) and recall that R^R = I, all the matrices R are cancelled, 
and we get 

G = 5: = tr (,„A') 

k>p 


where 


From (ITBl) we eventually have 


2k-l-N^ 

\do)kp — 9 ^kp 


J [X] = exp j J ti {rioX(t)) dt 

.0 


5 Isotropic processes 


Now we express the probability functional in the form 


F [A] = exp I - y L (A, dtA, d^A, ...) dt j , 


where L is the Lagrangian of the process. 

The change of variables ([9]) transforms F[A] DA to Px[X] DX where 

T T 

Px |A] = P J = exp I - / Lxdt + / tr (ri„X) dt 


(18) 


(19) 


Here 

Lx = L (F^XF, dt (F^XF) , (F^XF) , ...) 

contains generally not only derivatives but also integrals of X{t), since F[X] contains the 
T-exponent fllUp . Hereafter, we restrict our consideration with the isotropic processes: 

F [O^AO] = P [A] VO e SO{N) (20) 
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(This means that statistical properties of the process would not change under the global 
rotation of the reference frame.) We will show that for such processes, Lx does not contain 
the time integrals. 

The condition (12(11) implies that the matrix A and its time derivatives can contribute to 
L only in scalar combinations like 

tr {dt^A^)^^ ...) (21) 

But from ([8]) it follows 

dtR = OR , d'^R = {dtO + e^)R, dtR^ = R^O'^ ,...; 

accordingly, 

dt {R^XR) = R^ {e'^X + dtX + X9) R , d^{R^XR) = R'^ {...) R 

Hence, as we substitute R^XR for A in fl2TD . all the R^ are multiplied by R and vanish. 
As a result ) Lx contains scalar combinations of X and 6 and their derivatives, and does 
not contain R. According to (El), the matrix 6 is itself a function of X, so 

Lx[X] = Lx{X, dtX, d^X,...) 

Thus, the new Lagrangian Lx is a differential function of the matrix X. 


6 Isotropic (^-processes. 


There is an important particular case of absence of correlation between the values of A at 
different time moments: 

/ T \ 


P [A] = exp 


L (A) dt , 


( 22 ) 


\ 0 / 

the Lagrangian L being a function of A only, not of its derivatives. We also demand that 
A is isotropic (1201) . Then A can contribute to L only as a part of invariant combinations: 


L{A) = L (tr A, tr A?, tr AA'^, tr A^ 

The change of variables A i—)■ X results in the substitution X for A, since all the R and R^ 
vanish. Thus, 


Lx|A'] = L(21)U.^,. 


and 


tr {rjoX) dt j (23) 

To get the correlation fnnctions of X, one now jnst has to calculate the functional 
integral ffTTj) : 

{X,, (ti) ...Xkp (tn)) = N' J DX X,, (ti) ...Xkp (tn) Px [X] (24) 

where 

N' = 

Since X is an integration variable, it can be changed to A in the right-hand side of the 
expression. Note that the integral does not contain T-exponents. 



Px[X]= exp 


T T 

[ L{X)dt+ [ 


To simplify the calculations, one introduces the characteristic functional of a random 
process: 


Z[r]{t)] = { exp / tr {t] (t) A (t)) dt \)= DAP[A\exp / tr {r]A) dt (25) 


instead of P[A]. The averages can then be expressed as 


{Aij (ti) ... Akp (tn)) 


1 (5 

ZSrjij (ti) 


__ 6 __ 

^Vkp (^n) 


z 


r](t)=0 


(26) 


To determine the statistics of Q ([2]), (|1]), we have to calculate different correlators of X. 
From fl2T)) it then follows that for 5-processes the correlation function of X can also be 
written in terms of the same Z[ri\-. 


{Xij (ti) ...Xkp (tn)) 


1 5 

Z 5T]ij (ti) 


5 

^Vkp (dn) 


z 


v(t)=vo 


(27) 


The only difference from (|26|) is that the expression is calculated at the point ri{t) = rjo 
instead of rjit) = 0. 

It is convenient to introduce the generating functional for connected correlation functions 
dehned by 

Z[q] = (28) 
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The normalization requires hT[? 7 (t) = 0] = 0. Since Z[rj\ is a Fourier transform of -P[74], the 
isotropy (EO]) of P[A\ leads to isotropy of Z[rj\ and W[rf\\ 

W [0'^7]6\ = W [t]] VO e SO{N) 


Furthermore, if P[A] is a 5-process fl2^ . Z is also a continuous product of independent 
multipliers, and 

T 


^ [h (^)] = / w{rj{t))dt 


(29) 


The function wirf) is called a cumulant function (see, e.g., [H]). Via this function, one can 
calculate any of the correlators by consecutive differentiation: 

T 


{Aij (ti) ...Akp (tn)) — 


6 


^Tjij (tl) dTjkp (O 


exp 


w{ri)dt 


(30) 


ri{t)=0 


The first term of each correlator. 


d^w 

{Aij (ti) ...Akp {tn))c = -K - -K —(0)5(t2 - ti)6{t3 -h)... 6{tn - h) , 

OTJij . ••OT]}-p 

is called the connected correlation function and corresponds to the connected diagram (as 
ti,..., tn are represented by n points, and 5-functions make connections between them). 
The next terms contain the products of lower-order connected correlation functions (and 
hence, smaller sets of 5-functions) and correspond to non-connected diagrams. 

By analogy to fl271) . all the correlators of the X-variables can easily be obtained from 
the same expressions as these for A-variables by changing the point where the derivatives 
are taken. For the connected correlation functions we get: 


{Xij (ti) ...Xkp {tn))c — 


d^w 


dr]ij...dr]kj 


- ii)5(t3 - ^i) • • • 5(tn - ti) 


(31) 


7 The Lyapunov spectrum 

In this section we use the method described above to calculate the averages of the matrix 
elements X^j in the case of isotropic 5-processes. From (l3T|) it follows 

i^sq) = -^w (po) (32) 
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First, we show that the non-diagonal matrix elements are equal to zero. Actually, since 
w{r]) is isotropic, it depends on a combination of traces: 

w {r]) = w (tr T], trrf^ trr]r]'^,trr]^....'j (33) 

The derivative of each trace is 

^tr = Y, SisSn (h • • • 

Since the derivative is taken at rjkp = (r^o) kp ~ (see (CHD) , we get 6 isSjq6ij in each term, 
and non-diagonal elements vanish; 

d 

{Xsq) = ^- W {7]q) = 0 , s^q 

The averages of the diagonal components, 

Xs = {Xss) = (ps) {no summation) (34) 


are called the Lyapunov exponents. The set of is an important statistical characteristic 
of a process. It is used in many physical applications (e.g., to describe the separation of 
trajectories in a turbulent flow). 

Before we proceed to these calculations, we make one useful deduction: 

For the values A^, as well as other diagonal averages, the function w{ri) in ([3T]), fl3^ can 
be replaced by the diagonal cumulant function Wd where all the non-zero elements are set 
equal to zeros: 

Wd{vn,...,VNN) = w (35) 


This is caused by the diagonality of po so in fl24)) one can make an integration over the 
non-diagonal elements to get the ’diagonal’ probability function 

/ T 

^ ;dt 

k^p 


The cumulant function 


Pd [All, • • •, Ann] = / JI exp | - / T (A) i 
d k^p ' d 

corresponds to this PDF. 


7.1 Gaussian process 

Let the probability distribution of the process A{t) be Gaussian: 


P[A] = exp I / AijD.-j^^Akpdt 


(36) 
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where Dijkp is determined by the pair correlation function, 


2Dijkp5 (t - t') = {Aij (t) Akp {t')) 


Isotropy requires 


D'ljkp bSil^Sjp “ 1 “ cSipSjl^ 


Here a,b,c are constants (they must satisfy the condition P[A] > 0 for any A). The 
cumulant function corresponding to fl36|) is 


“(>)) = 'teAjtp’fc-p = a(trrif + btri)if + ctrrf 


(37) 


From (l3^ we then get 


= D (2s-I-N) , D = b + c 


(38) 


We note that the Gaussian spectrum is antisymmetric relative to the change s ^ iV —s + 1; 
in particular, 




0 


In many applications one needs the additional restriction trA = 0. From ((trAY) = 0 
it then follows ^ 


N Oj T 6 T c — 0 


Since a does not contribute to A^, this condition does not affect the spectrum (l38ll . 


To describe the statistics of Xgs = ps more accurately, one can also calculate its mean- 
square deviation from the average: denote 

Ps Xg , 

then 

(C. (G) iq (G)) = {ps (H) pq (h))^ = {ps (h) pq (G)) “ XgXg (39) 

According to (1^ . 

d d 

{is(tl)iq(t2)) = Dsq5(tl-t2) , ^*9=0 - O-^ (^o) (40) 

^Vss 

^The matrix then becomes singular, which formally corresponds to appearance of the multiplier 

Y\5{trA{t)). 

t 
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Simplifying the calculation by using only diagonal components of w{ri) in accordance with 
fl35|) . for the Gaussian function flTTI) we have 

= 2 {D5sq + a) 

Unlike the Lyapunov exponents, the dispersion of ps depends on the distribution of trA. 
In the case of traceless Gaussian matrices we get the well-known relation derived by na: 

(ef («i) ?? («2)> = w - Aj s («, - fa) 


7.2 Non-Gaussian process 


Let now A{t) be a non-Gaussian isotropic 5-process. As we have seen in the beginning of 
the section, for diagonal components of {Xgq) one can use fl3^ with w{ri) replaced by its 
diagonal part (|35]): 


A., = 


d 


drjs 


■Wd {r]ii... Vnn) 


(41) 


n=m 


The isotropy condition 
with different n: 


means that Wd must be a function of combinations 'Yh iVs 


Wd 


iVu • • • Vnn) = / I 5^ iVss ), ’ • • •) 

\ s s s / 


(42) 


Decomposing Wd into symmetric and antisymmetric (relative to rj ^ —rj) parts, we get 
(h) = {p) + {p) , (p) = ^ ^ 


One can see from fl30ll that the symmetric part of Wd contributes to the connected even 
order correlation funcitons of A, and wj contributes to those of odd orders. Furthermore, 
we decompose 

As = A^ + Aj" , (43) 


here A+ is produced by the antisymmetric part of Wd, and vice versa: 


A 


± 

S 



VO 


(44) 


Then, symmetric properties of Wd cause the symmetric properties of A* relative to po —po-, 
i.e., s ^ N -\-l — s\ 


A 


ib 

7V+1-S 


itA 


zb 

S 
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Hence, the symmetric part is determined by the odd-order correlators of A, and the 
asymmetric part depends on the even-order correlators]^ This is an important feature 
of the random processes: e.g., in theory of turbulence it is connected with the asymmetry 
of a turbulent flow relative to the change of time direction ([HI fin]). 

If the matrices A (and hence, X) are assumed to be traceless, we have 


s s s 

Generally, the Lyapunov spectrum is determined by the choice of / in fH2|) . Here we 
analyze one particular (important for physical applications) case of traceless 3x3 matrices 
A(t) with distribution close to the Gaussian. 

For N = 3 and trA = 0, the Lyapunov spectrum takes the form 


Ai = -A-f, 

A2 

A,s = A - 


(45) 


'3 - - 2 

According to (145]) . A 2 is determined by the odd, and A by the even part of Wd- If we restrict 
f|42l) to the hrst two terms, the traceless condition 


gives 


Wdiv) = T) ( ^ ^ r]qq 

<? \ q 


F{(Ev^J(Evqq)-Ev, 


Iqq 




(46) 


The coefficient F is called the asymmetry coefficient of the process; in accordance with 
(]30l) . it determines the third-order correlation of A, for example 


(All (H) All (^2) All (G)) — ^ {ti — ^2) S {ti — G) 


From (1411) it follows 

So, 


X, = 2D{s-2) + 2F{l- -{s - 2y 


A = 2D 


Ao — 2F 


(47) 


^Note that in the Gaussian case the Lyapunov spectrum is odd because there are no non-zero odd-order 
correlators. 
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For the second-order correlator ([39]), fHOj) we obtain 

The hrst term in this equation is Gaussian, the second term corresponds to the main non- 
Gaussian contribution. 

In non-Gaussian process, there is also the third-order connected correlator: 


Dsq = 2D I - 


'sq 


+ 2F {s + q - 4) i 1 - 


(6 ih) iq (G) ip (G)) = {ps (G) pq ih) Pp (G))c = F,qp S (tl - t 2 ) 5 {tl - G) 


where 


F = 


d d d 
dpss i^Vqq i^Vpp 


Wd ivo) 


Substituting fH6l) for w, we get 


Fsqp 2F ( 6sa T ^sv T ^ov 3 (5. 


'sq 


sp 


■'qp 


-'pq^sq 



{no summation) 


The hnite-polynomial approximation of Wd has a serious defect: such Wd corresponds 
to the probability density that is not positively dehned m)- But still it is a useful 
simplihcation, since the rest of the series does not make fundamental changes to lower 
-order correlators. 

There is a simple way to estimate the validity of cutting the third term in (06]). It is 
known that the Lyapunov indices must be ordered ([3l[T2]): 


Ai < A2 < A3 


From f05l) . f07|) we then have 

ir., 2T) 

\F\ < — 

' ' 3 

This is a necessary condition for cutting the series (06|) . In theory of turbulence, F is also 
required to be positive to provide the right sign of energy dissipation (which is associated 
with the third-order correlator of A). 


8 Statistics of the T-exponential 

In the previous section we discussed different moments of the matrix X. Now we proceed 
to the averages of its exponentials ([7|). As we mentioned above, the matrix R (which is 
the T-exponential of 6) separates from the other variables (IH|). The calculation of z (which 
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is, roughly speaking, the T-exponent of () is simplihed by the fact that = 0, so all the 
series are hnite. Here we analyze more accurately the statistics of d. According to (j7|), 

ds (t) = exp Ps (t) = exp ^ j Xss (r) dr 

Thus, its moments are 


((d. = ( exp I n y X,, (r) dr | ) (48) 

These characteristics are important in the applications; in turbulence, and in particular in 
the theory of passive scalar advection, they describe the separation of trajectories of liquid 
particles. 

Consider the characteristic functional and the cumulant function of the X{t) process: 

z'^-' !„(«)] = exp ^ J w^^\ri)d^ = ^exp ^ J tr {ri(t)X(t)) d^ ^ (49) 

According to fl22l) . ([23]), fl25l) they are related to those of the process A{t) by 

[r]{t)] = Z [r]{t) + T]o] Z~^ [r]o] , w^^\ri) = w{r] + r/o) - w(ho) 


In the previous section we saw that to calculate the moments of X, one can reduce w{r]) to 
the ’diagonal’ function Wd (?7ii, ■ ■ ■, ? 77 VAr) (13^ - Similarly, to calculate fl48|) . it is enough to 
consider 

wi^\Vs) = summation 

Then from (15^ . (140]) . fl4T]) it follows that can be expanded into the series: 

i.Ps) = Kvs + + ^Vs + ••• 

where the hrst coefficient is the Lyapunov index fl34p . the second is the covariation matrix, 
and the next coefficients correspond to higher correlators of X. 

From fl48l) . fl4^ we have 

{{ds it)T) = exp {wf\n)t) = exp (50) 

We see that the moments increase exponentially with time. From (150|) it also follows that 
in calculation of moments of the exponentials fj48l) . one can not neglect the contributions 
of higher-order connected correlators even if all the moments of X are dominated by non- 
connected diagrams, he., even if the integral J Xdt satishes the condition of the central 
limit theorem and acts as a Gaussian, exp J Xdt is still essentially non-Gaussian. 
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9 Smoothed variables 


The cumulant function determines completely the statistical properties of a 5-process. How¬ 
ever, it cannot be measured in experiments (though the statistical averages, which are its 
derivatives, are measurable), nor is the probability density functional a measurable value. 
(Even more, the PDF is, as opposed to w{r]), badly defined in the case of non-Gaussian 
5-process, it requires renormalization.) In this section, we introduce smoothed random 
variables which can be measured in experiments. We show that they appear to have same 
cumulant function as the 5-process, and their averages are related strongly to the correlators 
of A and X. 

Consider a random matrix 

1 



0 

The corresponding characteristic function is 


z {t]) = (exp (tr (^t]A))') (51) 

It is related to the characteristic functional Z [rjlt)] and its cumulant function w{r]) by 


z{t]) = Z [rjO (1 — f)] = exp / w {rjO (1 — t)) dt = exp w {t]) 


(62) 


So, w{7]) is the cumulant function for both the stochastic process A{t) and stochastic 
variable A. The averages 

can be obtained from the corresponding correlators fl30l) by omitting the delta-functions, 
the terms 

are equal to the coefficients of the connected correlation functions. 

One can express w{r]) by means of the probability density p{A) : from (|5T]l . (15^ one 
can easily derive 

w (rj) = \n J dAp (H) exp (tr (pA^) 

For example, if A{t) is a Gaussian process fl36l) then p{A) is Gaussian, too: 


p (A) = exp (^-^AijD.jlpAkp^ 
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and corresponds to the same cnmnlant fnnction (l37D . In the general case, p{A) can be 
written in the form 


V (^) = exp - V (A) 


where the fnnction V (A) corresponds to the non-Ganssian part. Then the cnmnlant fnnc¬ 
tion is 


w (r]) = In 



where exp v^jj nnderstood as a formal series. This allows, in particnlar, 

to restore the cnmnlant fnnction by known coefficients in connected correlators. Now as we 
have the smoothed variables, the cnmnlant fnnction is not necessary: p{A) contains all the 
information abont the delta-correlated random process A{t), and 


{Xij{ti)...Xkpitn))c — N 


dAAij...Akpp{A) exp {tr{poA)) 6{ti -tn) , 


In particnlar. 


N=( dAp{A)exp (tr{pQA)) 


-1 


Xs = N dAAssp{A) exp (tr{r]oA)) 


10 Conclusion 

In this paper we present a very simple technics to calcnlate the averages of time-ordered 
exponentials of random N x N matrices (E]) and, thns, describe the evolntion of linear 
stochastic systems. 

The simplihcation comes from changing the variables (|9]), (E]) that allows to separate the 
rotations of the eigenvectors and therefore to ’stabilize’ the rest of the solntion. The formal 
solntion ([2]) can then be rewritten as a fnnctional integral flTT]) with a simpler ’rotational’ 
T-exponent (ITnll . 

If the random process is isotropic, the ’rotational’ part can be excinded, and the av¬ 
erages (ITT]) become fnnctional integrals withont time-ordered prodncts. In the case of 6 - 
correlated in time process we calcnlate the Lyapnnov spectrnm fl3T|) and other correlation 
characteristics of the T-exponent fl3Tl) . 

For Ganssian probability distribntion of the matrices, the resnits coincide with those 
obtained earlier by [1]; the non-Ganssian case is analysed in details. We also describe 
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the relation of the statistics of the T-exponentials to the statistics of ’smoothed’ variables. 
They appear to be determined by same cumulant function, so that, measuring probability 
density function (or correlators) of the ’smoothed’ matrices, one can easily calculate, e.g. 
the Lyapunov spectrum. 

The important feature of the non-Gaussian probability distribution is that the higher- 
order connected correlation functions cannot be neglected when counting the exponential 
averages, even if the difference from the Gaussian is small and higher-order correlators 
of the exponent are dominated by the Gaussian (non-connected) contribution. One more 
difference from Gaussian is that the Lyapunov spectrum of a non-Gaussian process is, 
generally, asymmetric. This is very important for applications. In particular, in [13] the 
observed scaling space distribution of velocity and statistical properties of a turbulent flow 
were derived from the statistics of velocity deformation tensor at large scales based on a 
stochastic analog to Euler equation. 

It was shown that symmetry of the Lyapunov spectrum corresponds to time invariance 
of the flow, and hence must be broken in real flows. The right sign of energy flux (from larger 
to smaller scales in 3'’*) requires (trA^) < 0. Thus, one cannot restrict the consideration by 
Gaussian approximation: non-Gaussianity is of crucial importance. The results achieved in 
this paper allow to calculate the Lyapunov spectrum based on experimental measurements 
of large-scale velocity statistics, then derive the statistic characteristics of velocity held at 
small scales in accordance with [13], and compare them to the observations. 

This work is supported by the RAS program ’Nonlinear dynamics in mathematical and 
physical sciences’. 
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